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Using non-perturbative lattice method we studied hadronic screening correlators above and im- 
mediately below the deconfinement transition temperature, T c , in the quenched approximation with 
lattice spacing of 1/(8T) using clover improved Wilson fermions. Simulations were performed at 
temperatures T/T c — 0, 0.95 and 1.5. Mesonic screening correlators show no statistically signifi- 
cant thermal effects below T c , and clear evidence for weakly interacting quarks above T c . Baryon 
screening correlators yield similar physics above T c , but show precursor effects for chiral symmetry 
restoration below T c . 
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I. INTRODUCTION 



Experiments with relativistic heavy-ion collisions are beginning to probe the finer details of non-perturbative pre- 
dictions of QCD This is possible because fireballs produced in these collisions come close to thermal equilibrium 
0, and lattice QCD techniques can be used to predict expectation values of observables under these conditions Q. 
The question of why the fireballs come into equilibrium, i.e., why relaxation times are small remains outside the 
direct purview of lattice computations for now. Even so, one can ask why thermodynamics is applicable |5{. The 
answer can be found by testing whether any static correlation length in a QCD medium, £, (which is the inverse of a 
screening mass, fx) is much smaller than the size of the fireball. 

The only correlations which persist in the equilibrium thermodynamic system are static spatial correlations. These 
correlation functions are measured by introducing static probes into the equilibrium plasma and measuring the re- 
sponse of the medium. This response depends on the quantum numbers carried by the probes; so one can classify 
static correlators as glueball-like, meson-like and baryon-like probes with the usual quantum numbers of these quanti- 
ties corrected for the fact that the static spatial symmetries are different from the Poincare group. The group theory 
of these correlators was worked out for glueballs [6-8] and mesons [9( . In Appendix [X] we extend the representation 
theory to spin- 1/2 baryons. 

Meson- like screening masses have been studied in QCD in great detail [Iol - fl7l |. While we continue these studies 
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2.20 
1.99 


0.97 ±0.02 
0.94 ± 0.03 


1.05 ±0.05 
1.07 ±0.04 


0.86 ±0.16 
0.89 ±0.22 


0.75 ±0.10 
0.68 ±0.12 



TABLE I: Thermal shifts in masses and mass splittings at T = 0.95T C ; fin denotes the screening length in the hadronic channel 
H, and mn the mass at T — 0. We find that [in = ran within 95% confidence limits, indicating no significant thermal effects 
in the vector meson, V and nucleon, N+. Nor do we observe a thermal shift in the mass splittings between the V and the AV 
(axial vector) mesons. However, the splitting between the JV+ and its parity partner, N- , changes at finite temperature, and 
is a precursor to the restoration of chiral symmetry before the QCD phase transition. 
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here, the main focus of this work is the baryon sector. Baryon-like screening correlators were studied in the high- 
temperature phase dSEII) but there have been no follow-up studies in recent years, nor have these been extended to 
the low-temperature phase till now. Since these provide important inputs to the study of baryon-number fluctuations, 
which is of much interest in the experimental search for the critical point of QCD, we reboot the study of these 
correlators. 

Above T c we verify that baryonic screening lengths support the evidence for weakly interacting quarks which has 
been gathered over the years with meson screening correlators. Below, but very close to, T c we find no evidence of 
any finite temperature effects on the masses of the lowest mesonic resonances. However, we find that in the nucleon 
sector the mass of the opposite parity nucleon moves closer to the ground state (see Table HJ . These observations not 
only constrain models of quantum hadrodynamics [l8] , but could also have implications for the analysis of heavy-ion 
collision data. 

The paper is organized as follows. In Section |TT] we discuss technical details related to the simulation, inversion 
and analysis, including the details of interpolating operators used and the asymptotic fit forms. In Section IIIII and 
Section|lV]we display the results in detail. Section \V\ summarizes and concludes the work done in this paper. Various 
technical details arc presented in appendices. 

II. RUNS AND MEASUREMENTS 
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32 x 16 3 


6.03 


71 


1.7333 


0.1345, 0.1347 


0.95 


8 x 32 3 


6.03 


94 


1.7333 


0.1345, 0.1347 


1.5 


8 x 32 3 


6.332 


67 


1.5667 


0.1345, 0.1350, 0.1355 



TABLE II: The simulation and measurement parameters used in this work. /3 is the bare coupling for the Wilson gauge 
action, Csw is the clover coefficient (determined non-perturbatively in [3]), and k is the hopping parameter in the improved 
Wilson-Dirac operator. 

We measured correlation functions of operators with meson and baryon quantum numbers built out of clover- 
improved Wilson quarks [20l | in zero and finite temperature ensembles of configurations generated using the Wilson 
gauge action. The finite temperature configurations were generated on lattices with lattice spacing a = 1/(8T) with 
T = 0.95T C and 1.5T C . The spatial size of the lattice, L, was tuned such that LT = 4. We use the notation V = L 3 , 
V4 = V/T, N t = l/(aT) and N s = L/a. Zero temperature measurements were carried out on (2L) x L 3 lattices with 
coupling (3 = 6.03, which corresponds to the same lattice spacing as the finite temperature run below T c [2l|, [52]. We 
will describe later that this meant that L was 8.4 times the Compton wavelength of the lightest state on the lattice, 
as a result of which finite volume effects are under control. A zero temperature study at (3 — 6.332, corresponding 
to the thermal ensemble above T c was not made, but the corresponding choice of k was made so that we could use 
a previous study [23]. Details of the simulation and measurement parameters are tabulated in Table HU With the 
choices of k listed for /3 = 6.332, the values of m^/Tc are 1.52, 2.58, and 3.31 respectively. For the other parameters, 
the results of the mass measurements are described later. 

The gauge configurations were generated by a combination of one heat bath and three overrelaxation updates per 
Monte Carlo (MC) step. The first 2000 MC steps were discarded to ensure equilibration of the lattices. Subsequently, 
measurements were performed once in 200 MC steps, which is approximately 30 to 100 times the autocorrelation time 
of the action. As a result, it took 4 CPU hours to generate a single configuration on our largest lattice. 

The measurement of hadronic correlation function requires the inversion of the Dirac operator. This was done 
by a conjugate gradient algorithm which stopped when the residual dropped below (l/4/a 4 ) x 10~ 16 . Typically this 
was found to take about 2500 conjugate gradient steps, and about 6 hours of CPU time, for the smallest pseudo 
scalar mass used in our calculations below T c . The computations above T c are substantially cheaper; the inversions 
converged in approximately 0.67 CPU hour for our lightest pseudo scalar mass. 

The zero momentum correlation functions are 

S H (t)=J2( H H*,t)H(0,0)), where iJ(x, t) = $(x, t)T H ^(x, t), (1) 

X 

where ?/>(x, t) is the quark field at time t and spatial point x (we will also use the component notation x x = x, x 2 = y 
and X3 = z). The sum over all spatial sites of the point-to-point correlator projects on to zero spatial momentum. All 
Dirac, flavour, and colour indices are summed. T H is an appropriate Dirac, flavour matrix which gives the quantum 
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numbers of the meson H . In actual practice, we replaced the point source of eq. (JlJ by a wall source, in order to 
control the statistics better. This required gauge fixing; we fixed to the Coulomb gauge. The effect on the extraction 
of screening masses is discussed in Appendix [B] 

We used r ps = 75 for the isovector pseudoscalar (PS), T s — 1 for the isovector scalar (S), T v = 7$ for the vector (V) 
and r Av = 7^75 for the isovector axial vector (AV) channels. For the T — measurements, we considered propagation 
in the time direction as shown in eq. (JIJ, and since all the orthogonal directions are equivalent, we summed the V 
and AV propagators over the three polarizations i = 1,2,3. For the screening correlators all the polarizations of the 
V and AV are not equal, and we summed over only i = 1 and 2. 

The lattice mesonic currents need to be renormalized to connect them with the continuum currents. For quarks of 
mass m q , the lattice operators, i?(x, t), are multiplicatively renormalized to 

Z H (g 2 )(l + b H (g 2 )am q ^H(x,t), (2) 

where a is the lattice spacing, g 2 = 6/(3, and am q = (k — n c )/2. We use the results in the MS scheme where the 
renormalization constants are determined at the scale of 1/a. For non-perturbatively improved clover fermions, the 
factors Z V ^ AV and b v have been calculated non-perturbatively [2~jj . where the following interpolating formulae are 
found: 



1 - 0.7663g 2 + 0.0488g 4 _ 1 - 0.8496g 2 + 0.0610c/ 4 _ 1 - O.65I85 2 + 0.1226g 4 

1 - 0.6369 ff 2 ' Av ~ 1 - 0.7332 ff 2 ' v ~ 1 - 0.8467 5 2 ' 



For the other coefficients, we use the expression from the one-loop, tadpole improved perturbation theory [251 ] : 

Z PS:S = u (l - Zp S ,sg 2 ), (4) 



where uq is the tadpole factor, (P) = Uq, g 2 = g 2 ju%, and the one-loop coefficients z PS , s as obtained from [2(| are 
z PS = 0.107 - 0.019c sw + 0.0175| w and z s = 0.026 + 0.065c sw - - 0.012c| w where c sw = u%c sw . The one-loop order 
tadpole-improved forms of b H can be written as 

6 H = -(l + 6^ 2 ). (5) 

where one obtains from the calculations of [27j the values b l PS — 0.109, b\ = 0.070, and b\ v — 0.069, computed for 
the tree- level c sw = 1. 

Masses and screening masses can be obtained without renormalization. They were estimated, as usual, by fitting to 
an assumed cosine-hyperbolic form, and looking for agreement with effective masses. These imply assumptions about 
the spectral function at zero and finite temperatures, and we discuss these in more detail in Appendix [Bj Statistical 
errors on masses are obtained by a jackknifc procedure. The propagation of statistical errors is done by jackknife when 
they could be correlated, and by adding in quadrature when they are independent. The scalar propagator has strong 
quenching artifacts. This prevents the use of these correlators in estimating masses in the usual way (see Appendix 
0. 

Chiral symmetry restoration in the high temperature phase of QCD is signalled by pairwise equality of correlators 
which are related to each other by parity. The observation of some breaking of this symmetry has been an issue of 
some interest recently Here we introduce a new measure for this symmetry — 

" N t -1^ (S H+ (t) + S H _(t))> W 

where S H± (t) are parity partner correlators (for the screening correlators t is replaced by z and N t by N s ), and the 
angular brackets are averages over the gauge ensemble. The correlators at t = are left out of the sum to avoid 
problems with time doublers. We will use the convention of taking the positive parity partner for the label H. When 
chiral symmetry is broken we expect R H ~ O(l). 
For the nucleon operator we used 

JV a (x,t) = eabo{C^) P 6i>t(^tWp(^t)il> b e{x,t), (7) 

where a, /?, and 7 are Dirac indices, a, b, and c are colour indices, e is the Levi-Civita symbol, and C is the charge 
conjugation operator. The projection of the correlation function on to vanishing spatial momentum is performed 
by the usual means of summing over x at T = where we imposed periodic boundary conditions in all directions. 
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However, at finite temperature, anti-periodic boundary conditions must be imposed on quark fields in the Euclidean 
time direction. As a result, the screening correlator must be projected to the lowest Matsubara frequency [10]. The 
parity projection operators for the correlator whose propagation is measured in the direction fx is 

V£ = \{l±-f). (8) 

The group theory relevant for screening correlators is given in Appendix [AJ 

The zero momentum correlators of the two parities of the nucleon can be fitted to the behaviour expected of a 
Wilson fermion — 

S N+ (t) = c+ exp [-p N+ t] ± c_ exp [-p, N _ (N t - t)] 

S N _{t) = c-exp[-fi N _t] + c + exp[-p N+ (N t -t)}. (9) 

Because of the admixture from opposite parity states, these correlators change sign across the middle of the lattice. 
In this paper we used the absolute value of the correlators above. We note that R N does not require knowledge 
of Z N because the parity partners are generated by the same lattice operator. For the T = measurement, where 
£L N ± — m N ±. For the finite temperature measurement where we replace t by z, we have 

/4 ± =f N± - S in 2 (n/N t ), (10) 

when we use the prescription of [lOj. 

III. THE MESON SECTOR 





m n /T c m v /T c m AV /T c 


0.1345 
0.1347 


2.20 ± 0.03 3.55 ± 0.06 5.6 ± 0.4 
1.99 ± 0.03 3.53 ± 0.08 5.5 ± 0.5 



TABLE III: Meson masses in units of T c at T = and /3 = 6.03. 

The analysis of meson masses at T — is absolutely standard. The results are collected in Table InTl We choose to 
express all results in units of T c for two reasons. First, because quenched QCD is not open to experimental tests and 
hence quoting numbers in MeV units is based on assumptions which cannot be tested. We prefer to quote ratios of 
quantities which are computable in practice. Second, because in quenched QCD the critical coupling is known with 
high precision, and Nt = 8 is within the scaling region [281 ] , the statistical and systematic errors involved in using this 
as a scale are completely under control. 





H fi H /T c (J.H /m H 


2.20 


PS 2.18 ±0.02 0.99 ±0.02 
V 3.45 ±0.06 0.97 ±0.02 
AV 5.18 ± 0.09 0.93 ± 0.06 


1.99 


PS 1.97 ±0.02 0.99 ±0.02 
V 3.31 ±0.07 0.94 ±0.03 
AV 5.1 ±0.1 0.92 ±0.08 



TABLE IV: Meson screening masses, [i H at T = 0.95T C in units of T c and the corresponding T = meson mass, m H . 

The analysis of most mesonic correlators below T c is equally straightforward. The only subtlety has been mentioned 
earlier: since the zero momentum screening correlator is measured for separations along the z-direction, the three 
polarizations states of the V and AV are two spatial and one temporal. The temporal polarizations have behaviour 
distinct from the spatial polarizations @. We measure the screening masses of the spatial polarizations only. Our 
main results for the screening masses below T c are summarized in Table IIVI We discuss the relation between the 
screening mass and the pole mass in Appendix [Bj Our results show that the pole mass of the mesons is hardly affected 
by temperature. 
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m^/Tc 


T = 


T = 0.95T C 


Sint/m* X 2 /dof 


Eint/m^ X 2 / do f 


2.20 
1.99 


0.6 ±0.2 1.92 
0.5 ±0.2 1.28 


-0.1 ±0.8 0.25 
0.1 ±0.7 0.57 



TABLE V: -Ei n t in units of obtained by fitting the scalar correlator with the fit form in eq. (|C1|I . 



In Appendix [C] we discussed the reasons why scalar decays are visible in quenched QCD, and how the final-state 
interaction energy, i?i nt , can be extracted. The analysis involves fitting the functional form in eq. (IC1[) at both T = 
and for T < T c . The extracted physical parameters are given in Table El E- mt at T = is non-zero at the 95% 
significance level for both values of m v which we used. While the finite temperature measurements are statistically 
indistinguishable from these, they are consistent with zero within errors, due to the larger fit errors at T > 0. It is 
interesting to note that the central values of E^t at T > lie outside the 2a errors of the T = measurement. 

The SUi(2) x SUr(2) chiral symmetry of the QCD vacuum is restored above T c . The anomalous Uyi(l) remains 
broken asymptotically, although its effects are seen to be small at temperatures of 1.5T C . The most straightforward 
consequence of this is that the S/PS and the V/AV correlators should become pairwise degenerate in the high 
temperature phase. That this is the case is shown in Figure [U The figure also shows that the correlation functions 
are not too far from those expected in a theory of non-interacting quarks (also called the free field theory, FFT). 

At T = in the V/AV sector we found R v = 0.73 ± 0.02; at T = 0.95T C this drops marginally to R v = 0.59 ± 0.02 
for both the bare quark masses we have used. This indicates that chiral symmetry remains strongly broken up to 
0.95T C . A simple model of V/AV mixing below T c was presented in [l8j using a mixing parameter e, which can be 
adapted to our use by writing schematically 

S V (T) = (1 - e)S v (0) + eS AV (0), and S AV (T) = (1 - e)S AV (0) + eS v (0). (11) 

This gives R v (T) = (1 — 2e)R v (0). Using the values quoted above, we find e = 0.10 ± 0.02. In the unquenched theory 
one may expect much larger values of this parameter below T c [18]. We observed that at T = 1.5T C the value of R v 
drops as a power of the quark mass, vanishing as vanishes (see Figure [2]). However, there is no dependence of i?. s 
on the quark mass. 

We have extracted screening masses, fi H from the correlation functions. We found that effective masses showed a 
good plateau which agreed with fits, unlike previous experience with Wilson quarks at T > T c [TJ]. The reasons for 
this are explored in Appendix [B] where some effective mass plateaus are also displayed. Our results for the screening 
masses are collected in Table IVIII The pairwise degeneracy of the masses is quite evident. It is also evident that the 
screening masses are close to that expected in a theory of free quarks. The screening masses in the S/PS channel are 
7-9% smaller than that in the free theory, whereas in the V/AV channel they are within 2-3% of the free theory. 




FIG. 1: A signal of restored chiral symmetry above T c is that (a) the PS and S correlators and (b) the V and AV correlators 
become degenerate. The continuous curves are correlation functions in a FFT; in order to remove trivial artifacts, these have 
been computed on a lattice of the same size. 
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FIG. 2: Above T c R s is found to be almost independent of m„-/T c . However, R v goes to zero as a power of m w /T c , showing 
that correlation functions of the parity partners V and AV do become exactly degenerate in the limit. 



IV. THE BARYON SECTOR 



Analysis of baryon correlators requires more care than meson correlators. This is because states of opposite parity 
contribute to each parity projected correlator, as shown in Figure El As a result, both the fitting procedure and 
the extraction of effective masses is more complex than the analysis for mesons. In addition, the projection of the 
screening correlator on to zero momentum requires a twist to compensate for the thermal boundary condition. The 
subsequent extraction of a screening mass from the nucleon correlator also requires the subtraction of eq. (jlOjl , resulting 
in additional loss of precision. 

In Figure [3] it would appear that the finite temperature correlators are more nearly symmetric than those at T = 0. 
Whether or not there is an early onset of chiral symmetry restoration can be probed by constructing the measure 
R N , defined by eq. ([6]). At T = we find R N = 0.88 when m^/Tc = 2.20 and 0.89 at the lower quark mass. At finite 
temperature these change to R N = 0.8 and 0.83 respectively. Contrary to the visual impression created by Figure El 
once the covariances between the N± correlators are accounted for, the correlation functions themselves do not show 
any tendency towards early restoration of chiral symmetry. In fact if we adapt the model of eq. to this case, we 
find a mixing parameter e = 0.04 in the nucleon sector, which is even smaller than that found for the V/AV at the 
same temperature. 

The results of fitting the masses is shown in Table IVTl Here we indeed find some precursor effects of chiral symmetry 
restoration in the form of a thermal shift in the mass splitting between the baryon and its parity partner. We measured 
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FIG. 3: Nucleon correlators at (a) T — and (b) T = 0.95T C . Both sets are asymmetric, showing that the nucleon, iV+ and its 
parity partner, JV_ , are not degenerate. 
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m N+ /T c 


flN + /Tc 


m N _ /T c 


/ijv_ /T c 


2.20 
1.99 


5.10 ±0.10 
4.9 ±0.1 


5.3 ±0.2 
5.28 ±0.14 


6.8 ±0.4 
6.7 ±0.5 


6.6 ±0.4 
6.4 ±0.7 



TABLE VI: Nucleon masses at T = and screening masses at T = 0.95T C 
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FIG. 4: The analysis of mass splittings between iV + and N- for m^/T = 2.20 at (a) T = and (b) T = 0.95T C and for 
m^/T = 1.99 at (c) T = and (d) T = 0.95T C . The data points are effective masses, the band between the dashed lines is 
obtained by taking the difference of screening masses. The band between the full lines is obtained by direct fits to the splitting 
between the two masses, and gives the results shown in eq. (|12l) . 



the splitting by the ratio of the correlators, S N+ /S N _ , and extract the mass difference, Am, by fitting to an exponential 
form exp(— Ami). Since this method takes care of covariances between the correlators S N and S N _, we expect to 
control statistical errors better. The fit can be cross checked against the equivalent of effective masses for the ratio. 
Such checks are exhibited in Figured) The resulting values are 

m N _ -m N+ _ f 1,53 ±0.09 {m„/T c = 2.20) M»_ ~ M» + _ /l.l5±0.14 (m„/T c = 2.20) 

T c ~j 1.82 ±0.17 (m w /T c = 1.99) an T c ~j 1.24 ±0.20 (m w /T c = 1.99) ' ! '~ ! 

The thermal effect is significant, and leads to the results quoted in Table [IJ This is the most definite evidence to date 
about precursor effects to chiral symmetry restoration at T c . 

Above T c the correlators for N± become symmetric and degenerate. We find R N = 0.230±0.005 when m^/Tc = 3.31, 
R N = 0.146 ± 0.006 for m^/Tc = 2.58 and R N = 0.076 ± 0.009 for m„/T c = 1.52. These lead to a vanishing of R N 
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as a power of m w /T c . Chiral symmetry restoration is also seen in the fitted screening masses, displayed in Table \VU\ 
As in the meson sector, the screening masses are within 3-4% of those expected in a theory of free quarks. 



V. SUMMARY AND CONCLUSIONS 
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PS S V AV N + N_ 



FIG. 5: The ratio of screening masses measured at T — 1.5T C in quenched QCD with those in FFT. Although all the results 
are not taken for exactly the same quark mass, the effects of the quark mass are very small (less than 1% in this study). The 
results of these present computations are in agreement, within 95% confidence limits, with those from earlier computations 
with clover (13 ] and overlap [IB] quarks. However, the results from earlier computations with staggered quarks [TlHl3| differ 
considerably, being further away from the limit of free fermions. 



We have studied thermal correlations in quenched QCD in channels with various hadronic quantum numbers. The 
correlation functions are consistent with the usual picture of the QCD phase diagram — below T c long-distance corre- 
lations are mediated by hadrons and above T c these are mediated by weakly interacting quarks. It is approximately 
true that there are no thermal effects below T c (see Table Q} and that above T c the effects seen are consistent with 
almost non-interacting quarks (see Table EH]) . However, there are several details which may impact our understanding 
of the thermal behaviour of QCD. 

Above T c , at the temperature T = 1.5T C , we found strong signals of approximate chiral symmetry restoration in 
the near degeneracy of screening masses of hadronic parity partners, evident in Table I VIII In addition we found, as 
before [La - tla ]. that the correlation functions and screening masses are not far from those obtained in FFT, which is 
a model of non-interacting quarks. 

Since there is persistent confusion in the literature regarding the S/PS channels, we made a survey of the literature 
and found systematic effects which we summarize in Figure [SJ This work uses the smallest lattice spacing used in a 
quenched study of screening to date, and agrees with previous work using clover quarks 14j. It also agrees with results 
of studies using overlap quarks at smaller lattice spacing (l5j . However, it seems that studies with naive staggered 
quarks in the quenched approximation always indicate stronger deviations from FFT than any of these above studies. 



m-rr/Tc 


fips/T 




fly/T 




/Ujv + IT 


/T 


1.52 


5.54 ±0.02 


5.59 ±0.04 


5.88 ±0.02 


5.89 ±0.02 


8.72 ±0.10 


8.68 ±0.10 


2.58 


5.56 ±0.02 


5.59 ±0.02 


5.90 ±0.02 


5.91 ±0.02 


8.75 ±0.10 


8.73 ±0.10 


3.31 


5.61 ±0.02 


5.65 ±0.02 


5.95 ±0.02 


5.96 ±0.02 


8.82 ±0.08 


8.77 ±0.09 



TABLE VII: Hadron screening masses at T — 1.5T C for three different masses of quarks. In FFT all the mesonic (baryonic) 
screening masses are expected to be 5.99 (8.44), 5.995 (8.456) and 6.01 (8.48) on our lattices, when the bare quark mass is 
tuned to give rrw/Tc = 1.52, 2.58 and 3.31 respectively. 
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Since studies of screening with dynamical staggered quarks also claim such effects [12|, [16| , it is worth examining cutoff 
effects with staggered quarks more accurately in this context. 

Even if the screening masses are close to that in the free theory at T = 1.5T C , there are indications that there remain 
interesting physics in correlations. The deviation from FFT is the opposite of that seen in weak-coupling theory |29j. 
Using the measure R s , R v and R N of eq. ([6]), we examined correlations in parity partner channels as a function of the 
quark mass (parametrized through the T — value of m^). We found that as the quark mass is decreased R v and 
R N seem to vanish, as expected above T c . However, R s does not depend on the quark mass, and remains non-zero. 
This still does not resolve the question of the relative importance of the multiple possible origins of this effect [lj| , 
that will depend on future temperature scans. 

As highlighted in Table HI mesons and the nucleon ground state are indeed unaffected by thermal effects below T c . 
We also find that the splitting of the vector and axial vector meson masses is unaffected by raising the temperature 
to T = 0.95T C . However, the opposite parity excitation of the nucleon is seen to move closer to the ground state. A 
closing of the mass gap between these two would signal chiral symmetry restoration, and that is expected to occur 
only at T c . This is the first observation of precursors of the phase transition in the closing of the mass gaps of parity 
related states below T c . No such effects have been seen earlier either in the glue sector or with quarks, nor do we see 
any effect in the V/AV channels here. Since the parameters e introduced in eq. (|lf [) is much smaller in the quenched 
theory than expected in QCD with dynamical quarks, perhaps one can expect the change in the mass gap to be more 
significant in the full theory. 

A decreased splitting between the nucleon and its resonances implies that in a hadronic medium the interactions 
mediated through pions may keep the resonances in chemical equilibrium even after other baryons freeze out. This 
could result not only in changes in the net yield of nucleons, as seen at the LHC, but also possibly in the momentum 
distribution of nucleons. These implications clearly call for follow up lattice studies of nucleons below T c . 
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Appendix A: Splitting of Wilson-Dirac fermions at finite temperature 

The point group of a z-slice of the lattice at finite temperature is the group D4, which has 8 elements in 5 conjugacy 
classes. These classes are represented by the operators: the identity, n/2 rotations about the i-axis, R xy , 7r rotations 
about the £-axis, R xy , tt rotations about the y-axis, R xt , and ir rotations about the x+y-axis, R x+y = Rxy Rxt- Since 
the finite temperature theory is realized by periodic boundary conditions in the Euclidean time direction for bosons, 
the symmetry group for pure gauge operators is D\ = D4 x Z2(T) where T is the reflection operator in the xy plane 
Q. Since fermions are realized at finite temperature by anti-periodic boundary conditions, the Z 2 factor is not part 
of the symmetry group 

Acting on Dirac spinors, one can see immediately that 

R xy = e l ™ 12/4 = -^={1 +ia 12 ), R 2 xy =ia 12 , R 2 xt = ia 14 , (Af) 

where — [7^, 7^/2 and 7^ are the Euclidean Dirac gamma matrices. In the Euclidean chiral representation of the 
Dirac matrices, these three operators are block diagonal. The identity is always block diagonal, and R x+y is also block 
diagonal since both R xy and R xt are. As a result, under the action of D4, the 4-component Dirac spinor decomposes 
into two independent two-component objects, each corresponding to a specific helicity. 

The T — parity operator is P = (1 + 74)/2. Since this is (block) off-diagonal in the chiral basis, it mixes different 
helicity components. In the chiral symmetry restored phase, the two helicity components give degenerate correlation 
functions and one of the parity eigenstates vanishes. 



Appendix B: Screening Correlator and the spectral function 

Here we compile some relevant formulae for understanding the screening correlation function. We take the z 
direction as the direction of propagation, and examine correlators S(z) with z large. 

The screening correlator is easy to understand in terms of the transfer matrix of a z slice, and is therefore perfect 
for understanding the symmetries of the finite temperature transfer matrix. Information about thermal change of 
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properties of hadrons, on the other hand, is easier to obtain from the real-time retarded correlator j30j ] 

G R (x-T) = i6(x°){[H(x),H(0)]) T , (Bl) 

where Qx denotes thermal averaging at temperature T. The imaginary part of G R (p,T), the Fourier-transformed 
retarded correlator, is the spectral function a(p, T). A stable mesonic state, like the pion, contributes a term 

a H (p) ~ A H e( Po )S(p 2 ~m 2 H ). (B2) 

Here A H depends on the details of the operator. In particular, in our case of mesonic operators, A H depends on the 
nature of the source, and will be different for the wall and point sources. On the other hand, m H is the property of 
the hadronic state. 

At temperatures T sufficiently high but still below the transition temperature, one looks for thermal modifications 
of hadrons like a modification of the mass, or a thermal width. Such a modification may lead to, e.g., a relativistic 
Breit-Wigner form for the spectral function, 

a H ( P ;T) ~ Ze( P0 ) (|) (p2 _J^ +r2mr (B3) 

with both r and m possibly dependent on T. A similar form of the spectral function is introduced at zero temperature 
for a resonance. In that case this form is useful when the unitarity of the S-matrix is not an issue. When it is, then 
a more careful analysis is required where the decay products of the resonance are explicitly included. The treatment 
of the scalar meson in Appendix [Cl is an example. 

If we ignore such issues, then the screening correlator is connected to the spectral function, as usual, through the 
integral [3l| 

dp± rp. z r°° 2d po 



, 2tt J po 



It is easy to see that a change in the spectral function from the stable form eq. (|B2|) to a Breit-Wigner form eq. (|B3 
causes the large- z behavior of the screening correlation function to change to 



(B5) 



assuming T <C m H . One can see from this that the effect of a width as small as Y/m < 0.1, will be difficult to detect 
when the error bars are as small as a few percent. However, T/m ~ 0.2 should be visible, provided the spectral 
function of eq. (|B3[) can be used. 

On the other hand, the problem with the assumed Breit-Wigner form becomes clear if one examines effective 
masses instead of correlation functions. Then one finds that the effective masses would increase with z, contrary to 
the behaviour seen for the PS, V and AV. The reason is that in the lattice field theory, the unitarity of the S-matrix 
is enforced through the formulation. As a result, one cannot, in general, use a Breit-Wigner approximation to the 
spectral function in thermal QCD. 

Next we make a detailed comparison of the behavior of the screening correlator with the corresponding correlators 
in FFT. We used wall sources for FFT, and compared the effective masses, m(z), obtained on 32 2 x 8 x N z lattices 
(with N z = 80), with those obtained in the interacting theory. 

In Figure [S] we compare the PS and V screening masses for our lightest quarks (corresponding to m n /T = 1.52 at 
1.5T C ) to screening masses in FFT with similar bare quark mass. In the free theory, the screening correlator for wall 
source in these two channels are equal. The results for the FFT were obtained on a lattice with N z = 80 using the 
estimator lnC(z)/C(z + 1) for the screening masses. The plateaus in the effective masses shown is due to the use of 
a wall source. When a point source is used, neither the free theory nor the interacting theory shows a plateau 
These checks give us confidence that we are able to extract the asymptotic behavior of the screening correlator. 

The effective masses for the negative parity nucleon are also shown in Figure [5] In the interacting theory, a 
reasonable plateau is obtained in the effective mass, C(z) ~ exp(/i(z) x (N z — z)). The plateau is more pronounced 
in our data for heavier quarks. We see that the effective mass in the interacting theory is larger than that in the free 
theory. 

Since we use N s /N t — LT — 4 we also check finite volume effects in a theory of free quarks. In Figure [7] we show 
the screening masses for pseudoscalar meson and nucleon on lattices with N t = 8 and varying N s . With wall sources 
we find essentially no finite volume effect in either channels already for LT — 4. Since finite volume effects are larger 
in the free theory than in the interacting theory, therefore finite volume effects in our studies are expected to be small. 
In contrast, point sources give large finite volume effects [lij ]. 
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FIG. 6: Effective masses, in units of temperature, for meson and nucleon sources, at T — 1.5T C and in FFT. The black 
continuous curve denotes the FFT screening masses for Nt ~ 8. The black dashed line shows the continuum value, /i = 2nT 
[l^.l3^]. The data points and the bands show the effective masses and the fitted estimate for the asymptotic value, respectively. 
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FIG. 7: Mesonic and nucleonic screening masses in FFT as a function of Nt/N s ~ 1/LT for Nt = 8. The filled symbols are 
estimates from wall source quark propagators, while the unfilled symbols are estimated using point sources. In the latter case, 
since there is no plateau in the local mass, we use the prescription of [l4| and use the local mass at distance N z /4. 



Appendix C: Scalar Meson 




(a) (b) (c) 

FIG. 8: Quark line diagrams which couples with scalar propagator in the full theory. In the quenched theory (c) is absent. 

The non-trivial behaviour of quenched scalar correlators may be understood diagrammatically. In full QCD the 
scalar propagator has three kinds of topologies as shown in the figure above (33J. The third topology is absent in 
the quenched theory, as a result of which the scalar correlation function becomes negative when the quarks are light 
enough [33U35j ]. The missing diagram corresponds (for an isovector scalar) to the propagation of two different, but 
mass degenerate, pseudoscalars in the intermediate states; this is usually called a ttt]' (even in the two- flavour case). 
Due to the missing diagram, the quenched theory lacks reflection positivity. 

In full QCD one would be able to fit the scalar propagator with a simple cosh term. However, in the quenched 
theory the contribution of the missing term has to be cancelled (a quenched ghost) [33|. This leads to a complicated 
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spectral function resulting in the fit form which can be parametrized as [33- 35] 

S QxPT (t) = b{exp[-m s t} + exp[-m s (N t -t)}} 

-a s {(1 + rM) exp(-E n , w t) + (1 + m w (JV t - t)) exp(-E v ^(N t - t))} . (CI) 

The negative term is explicitly due to the quenched ghost. m n / s are the masses of the PS and S meson, E^^ is 
the energy of the ghost state (E V ' V = 2m n — E- U1 t) and — a s is the coupling to the ghost. Since the ghost exactly 
cancels a physical term, these two parameters are physical. For sufficiently small m w the second term dominates at 
intermediate distances, leading to negative values for the correlator. 

In our analysis of the scalar propagator we take the long-distance part of the propagator and extract from it both 
a s and E„i w . Using the renormalized propagator, we find the physical value of the coupling a s . These results are 
reported in Table [V] 

The spectral function of the finite temperature scalar screening correlator is not known. However, as for the other 
states, we will assume that the form of the spectral function for T < T c is the same as at T = 0, but the parameters 
could depend on T. With this assumption we will use the form in eq. ()Cip at finite temperature below T c . 
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